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Abstract 

With the use of theory developed earlier, bulk effects in ul- 
tracold neutron coherent inelastic scattering are considered 
both for solid and liquid target samples related to energy 
and momentum exchange with phonon and diffusion-like 
modes. For the neutron in a material trap, differential and 
integral probabilities for the energy transfer per bounce 
are presented in a simple analytic form which exhibits the 
parameter dependence. As an example, the theoretical val- 
ues for the ultracold neutron loss rate from a storage bot- 
tle with Fomblin coated walls and stainless steel walls are 
evaluated. Possible contribution from incoherent inelastic 
scattering on hydrogen contamination is discussed. 

PACS. 28.20.-v Neutron physics, 61.12.-q Neutron 
diffraction and scattering, 89.90.-|-n Other topics in areas 
of applied and interdisciplinary physics 



1 Introduction 

In the last few decades since the discovery ^ |2| of ul- 
tracold neutrons (UCN), great progress has been made in 
this field of physics. Now one considers UCN not only as a 
tool for studying fundamental properties of neutrons (life 
time and electric dipolc moment) but also as a method for 
investigating material surfaces and thin films E] . 

The main specific feature of UCN is their repulsion from 
material samples, which allows storing them for a long time 
in material vessels. This repulsion results from rescatter- 
ing of the neutron wave in the matter. Indeed, each atom 
in the target sample feels not only the incoming neutron 
wave but also the secondary waves from the other atoms. 
Effective secondary waves come mainly from surrounding 
volume ~ A^, where all ~ scatterers add coherently 
onto amplitude ~ n\^{b/X) (A is the neutron wave length, 
h is the neutron-nucleus scattering amplitude, and n is the 
density of scatterers). This amplitude may be neglected, as 
compared to the incoming neutron wave, if nb\^ ^ 1. This 
is just the cases for thermal and cold neutrons but not true 
for ultracold neutrons. Thus, for the ultracold neutrons 
rescattering in the matter is of crucial importance, and for 
the wave vector k < \J Airnb (if b > 0) the rescattering 
becomes the dominant process and results in the total re- 
flection from the surface. This dominant elastic scattering 



is simply described by a mirror potential U — 2TTh'^nb/m, 
where m is the neutron mass. 

Inelastic scattering of UCN on the trap walls also takes 
place and, in particular, (along with /3-decay and radia- 
tive capture) results in UCN losses from the material ves- 
sels. For a long time attention was focused on the inelas- 
tic transition of UCN to the thermal energy region (up- 
scattering), i.e. on large energy transfers (with probabil- 
ities per bounce 10~^ and higher). In the last few years 
the inelastic scattering with small energy transfers (" small 
heating" and " small cooling" ) was observed with probabil- 
ities per bounce lO-^-lQ-^ (see, e.g. El H El El ) . At the 
moment the data obtained for Fomblin oil as the sample 
appears to agree with theory. It was shown in 10 , the 
small heating influences substantially the UCN losses from 
the traps with walls coated by Fomblin oil. The data on 
solids like stainless steel, beryllium and copper are contro- 
versial, in particular, no small heating was observed in j^l 
for non-magnetic solid substances (stainless steel. Be, Cu) 
with the upper bound of ~ 10^* for the probability per 
bounce. 

Therefore, the inelastic UCN scattering is now of great 
interest, first due to the storage of UCN in material traps 
and, the second, for possible spectroscopic applications, 
e.g. for investigation of low frequency modes in surface 
layers and thin films 4 . On the other side, theoretical 
tools for the calculation of the inelastic scattering of UCN 
are rather undeveloped. One usually starts from Van Hove 
formula 11 for the differential cross section of neutron 
inelastic scattering on an ensemble of target nuclei. 



(1) 

Here and below we assume that k (e = h^k^ /2m) and k' 
(e' — h^k'^/2m) are the wave vectors (energies) for inci- 
dent and scattered neutrons, respectively, tkij = e — e' and 
Q = k — k' are the neutron energy and momentum trans- 
fers, R.i/(i) is the time dependent Heisenberg operator of 
the v-th nucleus position vector, and 6^ is the scattering 
amplitude on v-th nucleus. The averaging is understood 
over the initial state \i) of the nuclear ensemble. 

Note, however, that the Van Hove formula 1^ is based 
on the Born approximation (see, e.g. |12l 1130 and, hence, 
does not imply rescattering. Strictly speaking, it may be 
applied only to thermal and cold neutrons (sec also 14 ). 

Nevertheless, it is being used for UCN as the starting 
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point for some tricky ansatz. First, one notes that the 
inelastic cross section, both for coherent and incoherent 
scattering, predicted by is proportional to the total 
number of nuclei in the target. This allows to introduce 
the inelastic cross section for one nucleus. It should be 
stressed, that this definition is pure formal, because in fact 
the neutron is scattered by the whole sample as directly 
indicated by the presence of the correlation function in the 
right-hand side of 1^. Then the resulting cross section for 
neutron-nucleus scattering may be somehow extrapolated 
from cold to ultra-cold energy region (say assuming 1/w 
law). The sum of this inelastic cross section and that for 
the radiative capture on one nucleus is considered as the 
cross section of the total loss, which is supposed to define 
the imaginary part of the neutron-matter optical potential 
(see e.g. ^1 El)- Some model modification of the Van 
Hove theory for UCN was attempted in |17| . 

Note, that the observed UCN losses rule, 
higher than that theoretically expected. With this re- 
spect one uses the term "UCN anomalous losses" (see, e.g. 
|18l 1191 1^ ). In fact, while one speaks of a discrepancy 
between experiment and theory, there is really a lack of a 
reliable theoretical description for UCN losses. We would 
like to emphasize, that the rescattering results in a drastic 
change in elastic neutron scattering when one goes from 
the cold to the ultra cold region. Thus, one can expect an 
important infiuence of this transition on the inelastic scat- 
tering. Thus, a theory of the inelastic scattering of UCN, 
based on solid grounds is, certainly, needed. 

Recently, the authors presented a general theory of neu- 
tron scattering |14j . valid for the whole domain of slow 
neutrons from thermal to ultracold. The target sample 
is considered as a dynamical system with the eigenstates 
\j) and the eigenvalues ej. The theory is in fact a method 
of the solution of the -I- 1 body Schroedinger equation 
for a "neutron -I- target nuclei" system. The only approx- 
imation, which was used, is based on the fact that the 
neutron-nucleus potential is short-range (as compared to 
an interatomic distance and a neutron wave length) and 
deep (as compared to neutron and target nuclei energies). 

In this theory the amplitude of the neutron wave, 0, due 
to rescattering, is different at each nucleus, v, of the sample 
and depends on the energy already transferred, i.e. on the 
target sample state, j. In other words, the quantity (j>l is 
the neutron amplitude on the surface of the v-th nucleus 
provided that the target sample is at the state j. These 
amplitudes are determined by a set of linear equations. 
The differential cross section takes the form: 



k' 



dujdTt 27rk 

ly.iy J.J 

-too (2) 
X j e*'^*+*('^*~'^')*/''(j|e~**^^''(*)e**^^"'(°)|j'). 

— oo 

It was shown that for thermal and cold neutrons, when 
rescattering is not important, the equations for (j>l become 



very simple and result in i/)^ ~ Sijb^. Thus, the cross sec- 
tion takes the Van Hove form . For the opposite case of 
ultracold neutrons, when rescattering becomes the domi- 
nant process, we have solved the equations for the ampli- 
tudes (/)^. These amplitudes were presented as expansions 
in a small parameter Qu, where is the shift from the 
equilibrium position p^, of the v-th nucleus in the target 
sample (R^, = p^ + u^). 

The zero order in Qu corresponds to the scattering on 
the target sample with fixed, frozen (or infinitely heavy) 
nuclei and evidently can describe only elastic scattering. 
For inelastic scattering, one should consider the next orders 
in the small parameter Qu. It is convenient to introduce 
the renormalized neutron amplitudes = e^^P^(j)l/ f3^, 

and present them as series in Qu 



(3) 



Here (3„ is the scattering length on the bound nucleus [Tl] 
(this quantity slightly differs from 6^). 

Zero order amplitudes = ^{py) in the contin- 

uous media approximation were found to satisfy the 
Schroedinger equation with the optical potential U (p) = 
2Trh'^n{p)f3c/m, where /3c is the average ("coherent") scat- 
tering amplitude. Thus, the solution ipO^jP^) determines 
the zero order neutron amplitude on the surface of the i^-th 
nucleus inside a sample provided that the incident neutron 
has the wave vector k. 

Dynamical properties of the sample, and therefore the 
inelastic scattering, arise in the first order in Qu. The 
corresponding cross section was obtained in the form 
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(4) 

where ip{\i, v) = ■(/'(^k, i^) is to be interpreted as the wave 
function in the exit channel. Needless to say that the cross 
section Q should be averaged over target initial states \i) 
with thermal equilibrium density matrix. Below such av- 
eraging is implied. We emphasize that the approximation 
Qu ^ 1, which was used to derive Eq. is valid for a 
wide range of neutron energies, from the ultra cold to the 
thermal region. Thus, all transitions within this range can 
be investigated on the same base. 

In this paper we consider the case when the energies of 
initial neutrons are less than the barrier energy U . There- 
fore, the function ■(/'(k, p) decreases rapidly inward from 
the surface of the target, namely at the lenght on the scale 
of 10 nm. Then, clearly, the notions like the total number 
of nuclei-scatterers and the inelastic cross section for one 
nucleus are senseless. We are not introducing them. In this 
sense, our approach is quite different from that commonly 
used for UCN 
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We are not introducing and discussing the imaginary 
part of the optical potential, caused by the inelastic scat- 
tering. In our approach jj^, the functions of zero approx- 
imation, ipOi,!^) and ip{'k',v), are determined by optical 
potential with the imaginary part, defined only by radia- 
tive capture. At the same time, Eq.Q) gives the cross 
section of inelastic neutron scattering on the whole target, 
even on macroscopic one. 

Notice that in the paper we have performed straight- 
forward calculation of the total elastic cross section for a 
macroscopic plane sample and subbarier neutrons with the 
use of Eq.®. It was shown that it coincides with the 
cross section of the sample seen by incoming neutrons (as 
it should be). Thus, the probability of inelastic scattering 
per one bounce by macroscopic sample can be calculated 
as the ratio of the inelastic and elastic cross sections. Just 
in this way the corresponding probability is calculated in 
our approach. 

However, experiments with UCN in material traps do 
not allow a direct measurement of the differential cross sec- 
tion Q or differential probability per bounce. Parameters 
which can be measured are the total loss rate (in particu- 
lar, into the thermal region) and the energy distributions 
of pre- and post-storing neutrons. Since the probability 
of the small energy transfer is small, the multiple inelastic 
scattering of the same neutron may be neglected. There- 
fore, the energy spectrum of the stored neutrons can be 
directly linked to the total number of bounces with given 
material during the storage as well as to the probability 
of energy transfer per one bounce. In Section [21 we intro- 
duce the quantities measured in experiments, which will 
be analised and estimated. 

The purpose of this paper is to demonstrate a real pos- 
sibility to use our approach for analysis of the results of 
specific experiments on inelastic scattering of UCN. As the 
first step we assume that at the intrusion length of ~ 10 nm 
bulk properties of the matter are of dominant importance 
(discussion of the surface effects in the inelastic scatter- 
ing, e.g. contribution from visco-elastic waves for Fomblin 
oil or from hydrogen thin films, see in ^ 1^ I22|'l. 
Thus, in a uniform material the diagonal matrix element 
(il'u^t)u^, (0)|i) may exhibit a spatial dependence only as 
a function of — p^, and therefore allows a Fourier trans- 
form 



(^|<(^)ui,(0)|^) 



where Eq = / d^q/i^n f and J2u ^ I duj/27r. The corre- 
lation function il'-'Xq, w) for the target material is the most 
uncertain factor in our approach. We analize it in Section 

m 

Besides the usual phonon (sound) modes, we consider in 
details slow diffusion-like modes, namely, the contribution 
from thermo-diffusion and low frequency transverse modes 
in liquids. The idea is evident. Just the slow modes may 
be responsible for small energy transfer. However, for the 



sake of completeness, we perform calculation for all pos- 
sible final neutron energies, from zero up to the thermal 
region. Earlier, the role of the thermo-diffusion in transi- 
tion to the thermal region in solids was estimated in |17| 
and found insignificant. 

Really, for crystals, the dominating contribution to the 
up-scattering (into the thermal region) is due to phonon 
processes with the momentum transfer to the lattice. Our 
result for this contribution coincides with that obtained 
early in a model approach, quite different from ours (see 
e.g. 18,). 

All other processes, phonon- and thermo-diffusion ones, 
without lattice involvement, give small contribution into 
the scattering with large energy transfer. On the other 
hand, just these non-lattice processes are responsible for 
small energy transfers. In this paper, the corresponding 
probabilities for solids are calculated for the first time. 

Moreover, in liquids and amorphous matter, non-lattice 
processes are the only ones responsible for the energy trans- 
fers, both small and large. Thus, in this paper, the con- 
tribution of bulk modes both to small and large energy 
transfers for liquids and amorphous samples are estimated 
also for the first time. 

In the expression Q we take into account all terms with 
V ^ v' (as well as with v = ;/'), i.e. we calculate the coher- 
ent contribution to inelastic scattering. Thus, we use the 
average value (coherent scattering length) /3c instead 
in expressions for the inelastic cross section derived from 
Note, that to evaluate the incoherent contribution re- 
lated both to spin-flip and mixture of nuclei with different 
scattering lengths, one needs to separate additional term 
(proportional to /3|„c) from the general equation 10} keep- 
ing only V = v' . 

In this paper we restrict ourselves to the coherent in- 
elastic scattering. Thus, we realize that our consideration 
is unsuitable for treatment of neutron scattering on tar- 
gets with high hydrogen contamination. Our main goal is 
application to " hydrogen- free" materials. 

As the starting point we use the following expression for 
the coherent inelastic cross section fllj : 



da if, 



2mnk 



^i3"(q)i3^(q)f7'^(q,c.). 



(6) 



where 



B(q) = f3c 



V4V^(k',i/)V(k,z.)). (7) 



It describes the transition of an incident neutron with the 
wave vector k into an element d'^k' of the final wave vector 
space k'. This equation is valid for neutrons in a broad 
energy region, since the smallness of the parameter Qu, 
that was used for deriving Eq. , is valid even for thermal 
neutrons. Indeed, for the latter, the scattering on the op- 
tical potential (~ 100 neV) is not important. Then, if one 
replaces the functions ?/'(k, p) and V!(k',p) by the plane 
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waves e^^P and e^*"^ ^, respectively, then © transforms 
onto iQJ. 

In this paper we consider the inelastic scattering of initial 
ultracold neutrons to any final energies, below and above 
the potential barrier. Of course, to calculate the neutron 
inelastic cross section from ^ and {TJ one needs a specific 
model for the target sample to find solutions ^/'(k, p) and 
iplh'jp) for input and output channels. This problem is 
considered in Section^ 

In sectional it is shown that multi-dimensional integra- 
tion in @ and |(7J) can be performed explicitly with real- 
istic correlation functions. The details of the calculation 
are presented in Appendix. In section |S1 we analyze how 
analytical results for differential probability are changed 
when specified for the excitation mode and energy trans- 
fer. In section 13 numerical results are presented for two 
typical substances, stainless steel and Fomblin oil. Results 
obtained, both analytical and numerical, are discussed in 
section IHl Conclusion is given in section El 

2 UCN storage in material trap 

2.1 Probability for energy transfer per one 
bounce 

The probability for energy transfer per one bounce is nat- 
urally defined as a ratio of two cross sections, differential 
for energy transfer and the integral elastic one. The latter, 
in our case, is equal simply to the total sample area seen 
by the neutron, S± — Sk±/k — Scos9, 9 is the angle of 
incidence. With fixed initial wave vector k (or two its com- 
ponents fcj_ and fc||), the probability for the neutron after 
a bounce to have the kinetic energy e' is 



dw{k±, ki 



1 f da if. d^k' 



de' 



5, 



d^k' de' 



(8) 



In reality, the inelastic scattering of UCN with fixed 
k cannot be observed. One usually measures quantities 
somehow averaged over the initial momentum k. Let us 
introduce a quasi-classical distribution F{r, k) in position 
r and momentum k of UCN inside the trap, normalized by 
the condition 



where v — Kk/m is the neutron velocity, and mghs is the 
neutron potential energy in the earth's gravitational field. 
The total rate of inelastic transition to the energy interval 
dE' is of the form 

dN{E') 



— = ids / d^kv^Firs.liy, 

dE'dt J J 



k^>0 (11) 

d'w{k± , —fE' — mghs) 
^ dE' 

Let us assume that the momentum distribution is iso- 
tropic, i.e. i^(r, k) — F{r,k). Then, it is convenient to 
introduce a distribution /(r, e) in position and kinetic en- 
ergy of UCN inside the trap 



fir,s) = 27r[^y'v-sFir,k), k=^^, (12) 



normalized by the condition 



d'\ J def{r,e)^No. 



Then, the transition rate takes the form 



(13) 



dN[E') 
dE'dt 

where 



dS J devfirs,e) 





dw{e E' — mghs) 
dE' 



dw{e e') _ c?u)(fcj_, fc|[ — > e')^ 



de' 



de' 



(14) 



(15) 



is the differential probability of inelastic scattering aver- 
aged over the angle of incidence 9. 

Here and below averaging of any function f{kj_, k\\) over 
9 is defined as 

= ^ y d9sm9cos9 f{kcos9,ksm9). (16) 



/ d^kFir,k) = No, 



(9) 

where Nq is the total number of UCN in the trap. 

The number of neutrons, scattered by the element dS of 
the material sample at the position rs (and height hs) to 
the interval dE' of the total energy E' = e' + mghs during 
the time dt, is given by 

dN^ dE'dt Jd^kF{rs,k) (vdS) x 

dw{k±, fc|[ E' — mghs) 
dE' 



(10) 



In the simplest approximation we neglect the earth's 
gravitational field, and assume, first, the density nu of 
UCN in the trap is uniform over the volume and, second, 
the energy of UCN is fixed and equal eu, thus 



/(r,e) ^nuS{e- eu)- 
The inelastic transition rate is of the form 

dN(e') „ dw(eu e') 
= bnuvu- 



de'dt 



de' 



(17) 



(18) 



where S is the total area of the material surface inside the 
trap. 
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2.2 Probability for neutron losses from a 
trap 

When the final neutron energy E' exceeds the barrier en- 
ergy U, the neutron escapes from the trap. The probabihty 
of this event and the corresponding transition rate can be 
obtained as an integral over E' from U up to infinity from 
(|14|l (in the simplest approximation ~ as the similar inte- 
gral over e' from (llSfl ). 

The total escape probability is one of the quantities of 
practical interest. The second quantity of our interest is 
the small energy transfer. To be compared with measured 
values, both probabilities should be properly averaged over 
the energy and space distributions of the stored neutrons, 
and corrected for escape during the storage. Such averag- 
ing depends on the details of any specific experiment. On 
the other hand, the simple estimate (|18|l for the inelastic 
transition rate and its integral over e' allow useful rough 
comparison with measured values. 

3 Correlation function 

In this section we consider the correlation function fi*-' (q, ui), 
entering Eq.|(Bl). Correlation functions describe response 
(relaxation) of a substance after distortion of its statistical 
equilibrium by some external force characterized by lu and 
q. Their structures for solids and liquids are covered in 
many books and articles on fiuctuations and kinetics (see, 
e.g. ESllllESllini)- We restrict ourselves to the details 
necessary for what follows. 

Fast distortions (when ^ t^^, where r is some av- 
erage relaxation time) are relaxed by phonons (or sound 
waves). In this region the correlation function can be eas- 
ily obtained by expansion of the displacement vectors ul{t) 
in phonon amplitudes. The result is well known (see, e.g. 
Eq. (4.26) of Ref. |21|): 



(19) 



where the longitudinal and transverse parts are given by 

(20) 



and 



rii,t((j,cj) = 



pelt 



(21) 



(22) 



transverse (t) sound velocities Q.f and damping factors Ti^t 
to be different. 

For isotropic solids, the factors F^^^ can be taken propor- 
tional to the absorption coefficients of longitudinal and 
transverse sound (see, e.g. j^) 



cf p 



F 



1 
P' 



(24) 



(25) 



Here 7] and C, = Arj / ?i + Lp are the dynamic shear and lon- 
gitudinal viscosities {lp is the volume viscosity), v is the 
kinematic viscosity, a — CpjCy — 1, where Cp and Cy 
are the specific heats for constant pressure and volume, re- 
spectively, Ds — Si/ pCp is the thermo-diffusion coefficient 
(ae is the thermo-conductivity) . 

To describe a small energy exchange between UCN and a 
sample, the form of correlation function for small uj and q is 
of prime interest. Slow and long-range fluctuations are re- 
lated to coherent motion of a great number of media atoms 
(hydrodynamic modes). For such kinds of distortion, local 
statistical equilibrium comes first and macro-relaxation by 
hydrodynamic processes follows. Therefore, in the region 
of small UJ and q correlation functions for solids and liq- 
uids are governed by hydrodynamics and should depend 
on its parameters. Two different approaches, namely, phe- 
nomenological (from hydrodynamic fluctuations and Kubo 
theory [221 El) and quantum many- body theory _28_, give 
the same result. In hydrodynamics two kind of processes 
(and corresponding correlations) are distinguished, longi- 
tudinal and transverse. 

The longitudinal correlation function is related to " den- 
sity-density" fluctuations 



1 



2ttN 



( / ri(r' +r,t)n(r',0)dr') 



5'(q,tj), 



(26) 

where a Fourier transform, the "dynamical structure fac- 
tor", can be shown to be connected with r2;(q, w) by 



5(q,o.)c.^C!,(q',o.). 



(27) 



The vector q' in fi; is in the first Brillouin zone and q, 
the momentum transfer, may differ from q' by a reciprocal 
lattice vector. In the right-hand side of H27(l . we omit all 
terms proportional to 5{uj) and related to elastic scattering. 
Following |23[ I24| we get both for solids and liquids 



Here p — Mn is the material density (Af is the nuclear 
mass), and 



rio(w) 



T 



exp(?i|w|/r) — 1 h\uj\^T h\uj\ 



(23) 



is the occupation factor for a mode with the frequency uj 
at the temperature T. We allow the longitudinal (1) and 



Vti{q,uj) 



T 



1 



cjj.q^ - iujTi{q,uj] 



— c.c. 



(28) 



tpoj \UJ^ 

which looks like a phonon-type but with complicated form 
of "phonon absorption" 

r;(g,u;) 



p 



(29) 
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where Dt is the thermo-diffusion coefficient. The sub- 
scripts T and S indicate constant temperature and entropy, 
respectively, and 



4s 



= 1 



(30) 



For g2 < \u;\/Dt, Eqs. (jSHJ and result in the pure 
phonon-type correlation function with the sound velocity 
CIS and the absorption defined by (|24|1 . For arbitrary small 
CO and q (hydrodynamic region) 



1^1 g<-, 

7 7 

the function (|28|l takes the form 



(31) 



iT 

fli{q,uj) = — X 



Vc2j,(g2 + iuj/Ds)[q^ - uj^cf^ + i^J^l/cfs) J 

(32) 

It has two kind of poles, of phonon (p) and of thermo- 
diffusion (d) origin, and it is useful to separate them as 



niiq,u)^nP{q,uj) + nf{q,Lo). 
Here nf{q,uj) is given by with Tf (j^Hl and 



(33) 



nf{q,iu) 



2T 



pcfs\uj\ 



-p4 I 
^ d 



i ^ — 



^. (34) 



It is easy to see that additional diffusion-type pole in 
(|32|l originates from the pole in the "phonon absorption" 
Ti[q,w) (|^ . The complicated structure of F/ (9,0;) is in- 
trinsic only for longitudinal phonons. 

The transverse correlation function for solids has only 
one, phonon-like, pole. So, transverse excitations in solids 
are relaxed only by phonons, and the correlation function 
for small to and q has the same form (j22|l as for large q and 

OJ. 

In common liquid models, kinetics is assumed to go by 
finite jumps of atoms from one equilibrium position to the 
other after some average waiting time r. For high frequen- 
cies |a;| ^ a liquid behaves like a solid. Low frequency 
waves (1^1 ^ T^^) damp due to viscosity. This limit is re- 
produced by the correlation function for the liquid model 
with one relaxation time (see, e.g. j23j ) 



f7t(g,cj) 



T 



1 



where 



ipcjq^ \uj — iTt{q.,ijj) 



2 ? 

1 -I- iUT 



(35) 



(36) 



The relaxation time is usually estimated as r = rj/G — 
yjc^-, where r\ and v = rj/ p are dynamic and kinematic 



shear viscosities, respectively, G is the modulus of elastic- 
ity, and Cj = G/ p. 

The correlation function (|35|1 and H36|l exactly trans- 
forms to the form (|22(l with the damping factor 



M 

V 



Note, that for low frequencies 

|u;| « r-i = i 

V 

there exists a wide hydrodynamic region for q 



w Ct 
Ct V 



(37) 



(38) 



(39) 



where the Ea. H35|l (as well as H22(l ) transforms into the 
diffusion- type function 1)34(1 with damping factor (|37|l . It is 
of interest that this factor is proportional to v^^ in contrast 
to (E^ (see, e.g. ^). 

The correlation functions can be determined from first 
principles only for two limiting areas of variables, large 
and small w and q. In the intermediate area one should 
use some specific interpolation models for the substance 
considered. 

Thus, we use the correlation function as a sum of 
phonon- and diffusion-type terms. Contribution to the 
cross-section from each part can be calculated indepen- 
dently and summed afterwards. 

For the calculation of the cross-section we will need 
r2(q, Lj) as a function of {q^ = _|_ ^2^^ -g 
venient to present all cases considered above in the same 
form 



niqi 



pc^ 



(40) 



where 



n{qj 



(qi+p2)2 + r4' 



p 



^ii - (41) 



and a = e = 1 for all cases besides the longitudinal thermo- 
diffusion mode. In the latter case we have a = Cp j Cv — 1 
and e = 0. 

Note, that Fp (phonon absorption) is a small parameter 
for all values of uj and q that we are considering, and when 



0, 



%{q\_) 



(42) 



in agreement with the phonon correlation function usually 
used. 



4 Specific target sample: uniform 
thick layer 

In this section we will derive the function B(q) ((JJ) for a 
simple model commonly relevant to the UCN storage ex- 
periments, i.e. an uniform thick layer. The appropriate 
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choice of the model allows us to define the solutions ip(k, p) 
and '0(k', p) of the Schroedinger equation with optical po- 
tential for input and output channels. 

The problem is that the function v), being the neu- 
tron amplitude at the v-ih nucleus, is defined only inside 
the sample. But to get '(/'(k, v) in the continuous media ap- 
proximation, one has to find the scattering type solution 
of the Schroedinger equation with an external plane wave 
entering the sample, e^^P . For the input channel func- 
tion, ■(/'(k, i^), one may imagine the external plane wave 
as attributed to the incoming neutron. Let us assume it 
to come from the left. For the output channel function, 
"(/"(k'jp) = ?/'(— k',p), the corresponding external wave, 
g-zk is an auxiliary quantity and should be directed to 
the right, for back scattering, and to the left, for forward 
scattering. 

To include both cases, let us consider a uniform layer lo- 
cated at —a/2 < z < a/2 with the surface area S and the 
thickness a, much larger than the incident neutron wave 
length. The sample is supposed to have the density n and 
the scattering length /3c providing a constant optical poten- 
tial U . Note that in a general case, Pc is complex with an 
imaginary part that can be used to describe the radiative 
capture of UCN. 

We consider two solutions of the model scattering prob- 
lem, described above, with the left (L) and right (R) in- 
coming waves. Outside the target they have the form 



Re 



— iki 



z < -a/2, 
z > a/2. 



(43) 



, z < —a/2, 

jlf,ik±z^ z > a/2. 



(44) 

In both cases fcj^ = |k^| > is the normal component of 
the neutron momentum, and k|| is its component along the 
target surface. 

The functions inside the layer, determined by the bound- 
ary conditions, can be written (for constant U) as 



and the coefficients A are determined by 



The factor 



ika 



(48) 

7 = e'-" (49) 

oscillates very rapidly as a function of k± above the barrier 
(fc^ > fco); E^nd vanishes below the barrier (fc^ < fco). 

Below we use the "left" solution (|43|) for the input chan- 
nel function. Then, the output channels with the backward 
and forward scattering angles correspond to the "left", L, 
H43|l and "right", R, H44(l solutions, respectively. Inside the 
layer these solutions can be represented as 



X — - — y A. 

' I <j=±l 



1 



(50) 



All momenta and momentum dependent quantities H46|l - 
(|49() for the output channel are denoted by primes: k'j_, k' , 
t' , r' , A'^, 7'. They are interrelated by the same equations 
as that for the input channel. 

Combining (|45|l and H5()(l we obtain for Q 



p — i{k±+k'j^)a/2 



Here 



and 



Qii 



kl,. 



(51) 
(52) 

(53) 



where 



r/^ak + a'k', Er' = (77')'/'Ee^ 



(54) 

Summation over the nuclei on the plane parallel to the 
target surface gives 



1/2 
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(45) 



Sll = (27r)2n||^5(2)(Q||_q||_c;||), 



(55) 



where summation over the longitudinal reciprocal lattice 
Here k is the normal component of neutron momentum vector G|| arises for crystals, and n,, is the number of nuclei 



inside the target sample 



k=Jk{-k^, fco=47r/3cn. 



(46) 



per the unit area of the fixed plane (z = const). Therefore 
we have the equality 



Qii =qii +G 



(56) 



where ko is the value of the neutron wave number k±_ at the for the longitudinal transferred momentum. 



barrier, r and t are the refiection and transition coefficients 
for the target surface 



The cross section © is proportional to the second power 
of the quantity E||. Then, 



t = 



2k± 
kZ+k' 



(47) 



\^\= Si2TTfnlJ2s^'\Q\\ - - G\i), (57) 

Gil 
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where S is the surface area of the sample. 

For a non-crystalhne substance we replace summation 
over nuclear planes by integration 



a/2 



- (77') 



1/2, 



JiQi 



-a/2 



and obtain 



where 



= (77')'/VA(<7i-Q7'), 



A(x) 



2 sin(xa/2) 



(58) 



(59) 



(60) 



Here n± is the number of nuclear planes per the unit length 
along the z axis (nyn^ = n). The function A(x) peaks 
sharply for small value of the argument. 

For crystals the sum Yi"]^ has additional peaks when 
Q'\^ — q±^ is close to the transverse reciprocal lattice vector 
G\. In this case we take 



(77' 



Ml/2, 



Mqi- + G± - Qr ), 



(61) 



assuming the argument in A-function is in the first Bril- 
louin zone. Since the vector q±^ (originated from correla- 
tion function) is in the first Brillouin zone, the vector 
should compensate, if necessary, exceeding part of Q^'^ . 

5 Inelastic cross section and inelas- 
tic transition probability 

With the given correlation function H41() and quantity B 
(I51|l . the general expression for the inelastic cross section 
is fully defined and we can easily integrate over q|| with 
the help of (|57|l . Taking into account contributions from 
both, longitudinal and transverse parts of the correlation 
function, we find 



da'. 



Lt 



d^k' 



l\2 



(62) 



where the function 



l7l lYI 



|l-rVl'|l-^"7"l' 



(63) 



(T,cj' .r,r' 



includes all factors 7 and 7' that are strongly oscillating 
above the barrier and exponentially small below it. 

The factor J;,t(A, r\) depends on transverse momenta in- 
side the target in combinations 



yr 



(64) 



where A and 77 are assumed to be in the first Brillouin 
zone, and the difference from the full combinations 
and QY (see (|54|l ') is compensated by the transverse com- 
ponents G\ and G^ of the reciprocal lattice vectors. The 
factor J;,t(A, 77) is defined by the integral with the correla- 
tion function fli^t and has a different form for the longitu- 
dinal and transverse cases, 

4- 00 

Ji,t{X,v) = - / dq_LCli4q±)A{q_L - X)A{q^ - 'n)Fi,t{q±), 



where 

Fliqi 



(Qiiqii +Qr'gj-)(Qiiqii +Qr' q 



Ftiq±) = Ql + QT Q] 



Fl{q±). 



(65) 



(66) 



(67) 



Note, that the functions Fi,t{q±), as well as Ji^t, depend on 
indices ct, cr', r and r', but, to simplify notations, we omit 
them. Here and below we use q|| , defined from 



qii = Q|| - G||, 



(68) 



as a fixed quantity. The correlation function Q[q±), when 
parameterized as in (|41() by and p^, has a universal 
form for the longitudinal and transverse cases, and their 
specificity arises only from the numerical values of these 
parameters. Hence we shall omit the subscripts on fl{q±_) 
till the last, numerical stage. 

The inelastic transition probability (|HJl with the use of 
H62(l takes the form 



dwi^tik±,ku 



de' 



a(3ck no{uj) 



W^i,t(fc±,fc,| -.£'), (69) 



where the dimensionless factor W is given by 



n^k 



dk'^ / d(p|t'pn, 



i,t- 



(70) 



Here one integrates over k'_^ (up to k' = \/2me'/?i) and (p, 
the angle between k|| and k||, with fixed e'. 

For crystals, the final momentum k' is split by parts 
inside the first Brillouin zone, k'j^, and G. So the integral 
(O over k' should be understood as that over k'j^ and, if 
necessary, a sum over G. 

The detailed energy distribution inside a high Brillouin 
zone is of no interest, and we shall consider the full integrals 
over each zone, assuming that the final energy depends only 
on G, e' ~ h^G^/2m. In analogy with ©, ^ and ^ 
we define transition probability per bounce 



wi^t{k±,ku 



da'. 



i.t 



d^k' 



■d^k' 



(71) 
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which can be presented in the form 

wit{k±, fcii -> G) = X 

Mel, 

2m 



with fc'i 



where 



VK,,t(fci,A;|| ^G) = 



(73) 



is the dimensionless factor. 

The transition probabihties H69(l and H72|l are quite gen- 
eral and fully defined but rather complicated even for nu- 
merical analysis. However, a drastic simplification can be 
achieved by analytic calculations even for the general case. 

In our present publication we consider only a specific 
case that is when initial neutron energy is below the op- 
tical potential. Moreover, for the sake of simplicity in the 
main part of the paper we shall use only the results of the 
theoretical calculations referring to an Appendix, where 
the calculations are presented in details. Below we explain 
just the main steps of the solution. 

First, in A.l we calculate the integral H65|l. It turns out 
to be performed " almost exactly" by transforming (|65|l into 
a contour integral in the complex q± plane with the use of 
the four-pole structure of the correlation function (|41|l . 

The next step, calculation of the sum in H63(l . is per- 
formed in A. 2. Comparative importance of the 16 terms in 
the sum is very sensitive to the position of the initial and 
the final neutron energies. The result was obtained for the 
initial neutron energy below the barrier and for any final 
energy of the scattered neutron. After this stage, the in- 
elastic cross section is defined as a function of final neutron 
momentum, k'. 

The last stage, integration over angular distribution of 
the scattered neutrons, is done in A. 3 and A. 4. Here the 
momentum transfer between the neutron and crystal lat- 
tice is very specific and has required a quite different ap- 
proach (see A. 4). 

We write down the final results for the dimensionless 
quantities and (|75|l . 

For transitions inside the first Brillouin zone (G = 0) 
the factor l(7?Hl was found as 



Wi{k^,k\\ 



where 



Ak'Jkl 



k±f{k'j^)xo ^ 
37rfc 

2k^r^k'^ 

7rkaB{k'^-eLoycff^'' 
2k^f{k'^)xo 



Snk 



k'j_ < ka, 
A/{k'^ + ^k'l-^), k'^>ko, 



(74) 

(75) 
(76) 



fc'2 - kf., 



se 



^2 _ ^2 
"-0 "-J. J 



(77) 



(72) and 



M^) = \ n ■ i^^) 



The factors ((7^ for the case G 7^ were reduced to 

4fcj^a;o 



Wi{k^,ki\^G) = 
m(fci,fc|| -^G) 



3kee 
8k±xo 



(79) 



(80) 



The transition probabilities (O, ||7SJ) and lO, JHOI) can 
be further simplified by averaging over the neutron incident 
angle. There are two factors with angular dependence in 
(EH), IZgi and (Cnil, namely k±/{ka3) and k±f{k'j_)/k. Fol- 
lowing H16(l we introduce 



(81) 



where 



and 



1 



A{x) = (^arcsin y/x - x{l - x) ^ , (82) 



' k '^k^^U'W' 



(83) 



where the function B{x,x') has a different form in three 
regions of the variables: a;'<l,l<a;'<l + x, l + a;<a;', 
and can be presented as 



with 



B{x, x') = C(x, -) - Cix, (84) 

X X 



min(l,j/) 



Cix,y)=x J ^il-z)iy~z) dz. (85) 

miii(0,y) 

Note the two simple limits 

B{x,x') :^ x/2, for |x' - x| < 1, 

(86) 



B{x,x') ~ {l/S)./^^, for x'>l. 

The final expressions for Wi^t, defined in H7UI) and (|73|l 
averaged over the angle of incidence, are of the form 



xq £ e 



T^k" 



■ 2^k{k'^ -euj^/c^f^ "U 



(87) 



(88) 



The second terms both in (|74|l and (|87(l are valid only for 
large energy transfer (see H199f) , H2UU|) and text before) . 
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6 Physical application. Specific where v = \/2e/m is the initial neutron velocity, and 
ultracold neutron interaction with 
solids and liquids 



ci(y) 



1 V^c 



2.315. 



(96) 



In the previous section, specific cases of general equations 
were analyzed with the use of mathematical criteria. In 
this section, we consider the problems from a physical 
point of view, making classification by the relaxation pro- 
cesses, phonon-like or diffusion-like types, and the range 
of the energy transfer, with a special attention to the re- 
gion of small heating and cooling, as well as to upscat- 
tering region. The substances under investigation are of 
crystalline, amorphous and liquid type, commonly used in 
experiments. 

There are two longitudinal modes, phonon-like and 
thermo-diffusion, both for solids and liquids, and only 
one transverse mode. The latter is of the pure phonon 
type for solids, and of the combined type for liquids, that 
transforms from the diffusion type at small energy trans- 
fer to that of phonon type at large energy transfer (see 
Section OJ. 

6.1 Phonon-like mode, G = 

Phonon-like modes involve low damping, i.e. Tp —^ 0. 
Therefore, one can neglect the last term in H87(l and use 
for f7^ X(){i^) — |w|/c. For sohds, when G = 0, the sum 
over longitudinal and transverse modes gives from H87(l and 

(EH 



dwp{e e') _ (3c\uj\no{uj) ^ £ 

de' ~ ttMc^ ^U'u'' 



where c is defined by 



(90) 



(91) 



For the low energy part of the spectrum one can use (|23|l 
and find 

dl' -2^^^U'U^' k-e|«r. (92) 

This part of the spectrum is completely determined by the 
function B{x,x'). Here and below we use the quantities 



\/2toT 



Vt 



(93) 



The contribution of inelastic processes with G = to 
the total probability of the UCN escape from the trap can 
be estimated by integration of (|90|l from the barrier U up 
to the Debye energy 



Ed = hc{6TT^n 



,1/3 



Since U <^ en, it gives 



(95) 



For liquids, the contribution of the longitudinal phonon- 
hke mode is of the form (|^ and where 1/c^ 

should be replaced by l/3cf . The transverse mode for liq- 
uids is considered below (see Sec. 16.411 . 

6.2 Phonon upscattering, G 7^ 

The main process for upscattering in crystals goes with 
momentum transfer from the lattice. Assuming Fp — > 
and summing over longitudinal and transverse phonons, 
we obtain for the total upscattering probability from H72|l 
and Ip^ 



Am?Mc 



G||,Gi>0 



(97) 



where c is given by 1)91(1 . Restriction in the sum with > 
is made since the summation over backward and forward 
scattering was performed earlier. 

To estimate the probability we replace summation over 
> and G|| by integration 



E 



G|i,Gi>0 



d^G 



(2^/ao) 



(98) 



where ao = {jMlp is the lattice constant. Integral over 
e' = h^G^/2m up to Debye energy gives 

3 



V(£)=c,(-)^f^ 



where 



T 



C2(y) 



x^/'^dx 



1 y^c 



3.745. 



(100) 



The high energy upscattering due to phonons at G 7^ 0, 
was considered earlier ^1 by other approaches. The 
probability coincides with the result, presented, for 
example, in review paper |18| . 

6.3 Longitudinal thermo-diffusion mode, 

G = 

The differential probability of inelastic scattering via this 
mode, which follows from H69|) and H87|) with e = 0, can be 
presented in the form 



(94) dwd{e^e') aPcTdnoiuj 



de' 3^/2 7^Mcf 

SUTd e 
2Vme' U 




(101) 
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where the last term in (|87l) . vahd only for large energy 
transfer (e' ^ e), is written with factors 6'(e' — e) and (1 — 
g-e /lOe)^ which provide attenuation at small e' {0{x) — 
for X < and d{x) = 1 for x > 0). The damping parameter, 
r^, given in ll^ . can be presented as 



k 



|£ -g| 
due 



1/2 



where 



do = 



2mDs 



(102) 



(103) 



is the dimensionless thermo-diffusion coefficient. 

For small energy transfer, e' ^ e, the dominating first 
term in H101() takes the form 



dwd{s — > e') aPckj 



VVt 



de' 



QirMcf vl^/2dL 



(104) 



where vq = ^Jwjm is the barrier velocity. It is clear, 
that the divergence at e' ^ e is integrable. This result for 
longitudinal thermo-diffusion and for fixed components k±_ 
and A:|| of the initial neutron momentum was obtained in 

m 

For large energy transfer, the main contribution comes 
from the second term in H1Q1|I . It should be remembered 
that diffusion modes appear in hydrodynamic approxima- 
tion. So, the energy transfer via this mode has sense only 
inside the hydrodynamic region H31|) . where \uj\ <C cf/Ds- 
It gives the boundary for e' , where (|101|l is reasonable. 



, 2mc? 

do 



(105) 



If Ed < £d, the integration can be performed up to the 
Debye energy. It gives 



Wd{e) = ci( — ) 



T ' Ando M\ci I " U 



(106) 



6.4 Transverse mode for liquids, G = 

Taking a — e = 1 (see text after (|41|) '). we get from (|69ll 
and H88|l the differential probability of inelastic scattering 

dwt{e-^£') 2l3cXo{uj)no{uj) e e' 
= S.Afc? ^^C/'C/)' 

where XQ{bj) is given by H78() with c — > q, while the damp- 
ing factor Fj is of the form (|S7|) . 

Hydrodynamic region (j38(l . |a;| ^ c^/;^, corresponds to 
the final neutron energy 



, 2mc? 

di, 



where 



di, = 



2mv 



(108) 
(109) 



is the dimensionless viscosity. On the other hand, H1U8|) 
means F^ = \ijj\/v 3> /cl, which simplifies xt){uj), 



In the region e' ^ e, (|107|l takes the form 

dwt{e e') Pckr 
de' 




(110) 



(111) 



similar to H1Q4() for longitudinal thermo-diffusion. 

Beyond the hydrodynamic region, i.e. for large energy 
transfer, we have 



rf « 



Ct 



(112) 



and (|107|l transforms into (|90|l for transverse phonon-like 
mode, where should be replaced by 2/3C(. 

6.5 Relative contribution of phonons and 
thermo-diffusion to large energy trans- 
fer 

Phonon contribution to large energy transfer essentially 
depends on momentum exchange with the lattice. The 
probability for upscattering with (G 7^ 0, (|99|) ') and with- 
out (G = 0, (EHl) such exchange differ by the factor v/c 
(^ lO^'^ for UCN). The physical reason is simple. Inelastic 
scattering can be visualized as an absorption by the neu- 
tron of a quantum from the media with the energy huj and 
momentum hq. The requirement for energy and momen- 
tum conservation does not allow this process with "free" 
quanta. But for UCN the normal component of momen- 
tum is pure imaginary, zae, and instead of a (5-function, as 
in l|55|l for q||, one has from (|58|I - H6()|I for q± 



^i{q±+G±-(j'k')a/2 

sc + i{q± + G±- cr'k') ' 



(113) 



Integrals over q± with (|113|l are very sensitive to the 
minimum value of imaginary factor, q± + G± — a'k' , which 
may be reached in the allowed q±^ range. For sound-like 
excitations, a rough estimate gives 



k' <t:k', 

c c 



(114) 



where v' and k' are the velocity and wave number of scat- 
tered neutron. It is evident, that the large factor cr'fc' in 
the denominator of (|113|l can be compensated only by G± . 
Then H113|) can result in the large factor 1/ae, which is 
seen in (|7^ . l(5n|l . With G± = one may expect, instead 
of 1/ae, the factor 1/fc^, the result at least by the factor 
v/c smaller. 

Note, that the factor l/ai is of the scale of the UCN in- 
trusion length into the wall of the trap. Therefore, it is 
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natural that the inelastic scattering probability is propor- 
tional to 1/as. Thus, the factor l/k'j_ instead of 1/se points 
to substantial suppression of the large energy transfer by 
the phonon mode with G = 0. 

The thermo-diffusion mode can be also effective for the 
large energy transfer, but the physics is quite different. 
Parameters uj and q of " diffusion quanta" are not strongly 
coupled, as in the phonon case, (|114|l . Even with co fixed 
by the energy transfer, qj_ is still allowed to be varied up 
to ~ fc'. Thus, H113(l can reach 1/ae even for G = 0. The 
result can be seen from the last term in H74|l . It vanishes 
in the limit F — > (for phonons), but is proportional to 
1/ie for thermo-diffusion and results in the probability for 
upscattering H106() which may be comparable to that for 
phonons with G 7^ 0, 

7 Results for two specific target 
samples 

In previous section we have obtained analytic results for 
each type of relaxation modes. To get the physical result 
for a specific experiment, one should combine several an- 
alytic formulae with the relaxation modes and parameters 
appropriately chosen. To demonstrate this procedure and 
examine general trends of final results and their sensitivity 
to physical parameters, it is worthwhile to consider some 
specific examples. 

We have chosen stainless steel (SS) and Fomblin fluid 
(FF), two materials which are often used in real exper- 
iments and, on the other hand, have quite different re- 
laxation properties. Both materials are not of a simple 
structure, and we shall use for them simplified descriptions, 
which are specified below. 

7.1 Parameters for stainless steel 

The nuclei ^^Fe are the dominant ones in stainless steel. 
Thus, we take A = 56, where A = M/m is the mass num- 
ber of the scatterer. Then, p = 7.8 g/cm^ and, therefore, 
n = 8.4 • 10^^ cm"3 and ao = 2.3 • 10"^ cm. For the coher- 
ent scattering length we take be ~ 0.87 • 10^^^ cm, which 
gives for the barrier energy U — 192 neV. Sound veloc- 
ities at T = 293 K are equal to: c; — 5850 m/sec and 
Ct = 3230 m/sec. 

For the case of stainless steel one should consider con- 
tributions from longitudinal and transverse phonons and 
thermo-diffusion, the main part of the diffusion-like pro- 
cesses. For the latter we need a — Cp/Cy — 1 and the 
thermo-diffusion coefficient Ds, 



TK(3 



ae 



(115) 



where K = —V{dP/dV) is the volume elasticity, and 
Pp = [dV / dT)p /V is the thermal coefficient of volume ex- 
pansion. The corresponding values at normal temperature 



T — 293 K are listed in Table 1 as well as the calculated 
quantities a, Ds, and dimcnsionless parameters do and 

7.2 Parameters for Fomblin fluid 

Fomblin is a hydrogen-free fluorinated oil (it is known also 
as perfiuropolyether). It is quite viscous at room temper- 
ature and used to cover the walls of UCN storage traps 
because of the small neutron absorption and upscattering 
rate. There are different modifications of Fomblin oil and 
not all of them are described in the literature in detail. 
Thus we use typical values for Fomblin parameters. In ad- 
dition, we ignore complex structure of Fomblin molecules 
and are treating it as a sim pie liquid. 

The stoichiometry is roughly CsFgO. Thus we assume 
that the ^^F nucleus dominates in the fluid and A — IQ is 
the mass number of scatterer. We take p — 1.9 g/cm'^, then 
n = 7.15 • 10^^ cm^^. For the coherent scattering length 
we use he = 0.565 • lO"^'^ cm, which gives U = 106 neV 
in accordance with experiment. For all temperatures we 
take flxed sound velocities: q ~ 1900 m/sec and Ct = 
1500 m/sec. 

Parameters K, f3p, Cp, ae, v as well as the calculated a, 
Ds, and dimcnsionless quantities do and di, are presented 
in Table 1 for two temperatures: T — 293 K and 373 K. 
Among these parameters viscosity v is the most sensitive 
to the temperature. 

7.3 Numerical results 

We have calculated the probability of inelastic scattering 
for UCN with initial energy near half of the barrier height, 
that is for e = 80 neV for stainless steel and e = 40 neV 
for Fomblin. Contributions of phonons and diffusion modes 
are shown separately. Results for the differential proba- 
bility are presented on Figs. QEl and El 
Integral probabilities 



dw{e 



de' 



■de' 



(116) 



of transition to the interval U < e' < Sf per bounce are 
shown on Figs. ^ Eland [HI For Fomblin both differential 
and integral inelastic scattering probabilities are presented 
for two different temperatures: 293 K and 373 K. 

Test calculations have been carried out starting from the 
exact expressions H69|l and H70|l . Comparison with results 
from simplifled analytic formulae has shown quite small 
deviations in the whole energy range. So, we present nu- 
merical results based on analytic formulae. 

For stainless steel, the phonon contribution at G = 
was calculated from ()90|l and thermo-diffusion contribution 
- from (|101|l . Results at T = 293 K are shown by solid and 
dashed lines, respectively, on Figs. Q and 01 The flnal neu- 
tron energy for stainless steel is restricted from above by 
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Table 1: Parameters for stainless steel (SS) and Fomblin fluid (FF). 





T, K 


K, Pa 


dp, K-i 


g-K 


W 
m-K 


cm^ 
sec 


a 


sec 


do 




ss 


293 


170-10^ 


3.6-10-5 


0.45 


80 


0.01 


0.018 


2.28-10-1 


717.5 


31.5 


FF 


293 


3-10'' 


6-10-^ 
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Figure 1: Inelastic scattering for UCN of 80 neV on 
stainless steel at T = 293 K: the differential probability 
dwie/dsf per bounce as function of the final neutron en- 
ergy Sf. Solid line - phonon contribution, dashed line - 
thermo-diffusion contribution. 



Figure 2: Inelastic scattering for UCN of 40 neV on 
Fomblin oil at T = 293 K: the differential probability 
dwie/de f per bounce as function of the final neutron energy 
£/. Solid line - longitudinal sound contribution, dotted 
line - transverse sound contribution, dashed line - thermo- 
diffusion contribution. 
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Figure 3: Inelastic scattering for UCN of 40 ncV on 
Fomblin oil at T = 373 K: the differential probability 
dwie/dsf per bounce as function of the final neutron energy 
£f. Solid line - longitudinal sound contribution, dotted 
line - transverse sound contribution, dashed line - thermo- 
diffusion contribution. 



Figure 5: Inelastic scattering for UCN of 40 neV on 
Fomblin oil at T = 293 K: the integral probability Wie 
per bounce of transition to the interval U < e < Sf as 
function of £/. Solid line - longitudinal sound contribu- 
tion, dotted line transverse sound contribution, dashed 
line - thermo-diffusion contribution. 




Figure 4: Inelastic scattering for UCN of 80 neV on stain- 
less steel at T = 293 K: the integral probability Wie per 
bounce of transition to the interval U < £ < £/ as func- 
tion of £/. Solid line - phonon contribution, dashed line - 
thermo-diffusion contribution. 



Figure 6: Inelastic scattering for UCN of 40 neV on 
Fomblin oil at T = 373 K: the integral probability Wie 
per one bounce of transition to the interval U < e < Sf 
as function of £/. Solid line - longitudinal sound contribu- 
tion, dotted line - transverse sound contribution, dashed 
line - thermo-diffusion contribution. 
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the Dcbye energy Ed ~ 0.040 cV (I190|l with sound velocity 
c from H91|l . The thermo-difFusion contribution should be 
restricted from above by the limit energy Sd — 0.001 eV 
(see (pU). 

Integral probabilities on Fig. 0] calculated for G = 0, 
should be compared with the total upscattering probability 
at G 7^ 0, which is given by (|^ . At the temperature 
T = 293 K, it is equal to 



0.7-10- 



(117) 



For Fomblin, the longitudinal sound (phonon) contribu- 
tion is given by H90|l with l/3cf (see the end of 
Sec. I6.1|l , the transverse sound contribution - by H107|l and 
thermo-diffusion contribution - by l|101l) . They were calcu- 
lated both at T = 293 K and T = 373 K and are shown by 
solid, dotted and dashed lines, respectively, on Figs.|21|21ISl 
andini The final neutron energy is restricted from above by 
the Debye energy ed — 0.017 eV (|190|l with sound velocity 
from (PIJ. For Fombhn, the energy Ed = 0.023 eV (|105|) . 
the limit for thermo-diffusion contribution, is greater than 
the Debye energy ed- 



8 Discussion 

Differential spectra (Figs. ^ HI and O have well devel- 
oped maxima for both modes, phonon (with G = 0) and 
diffusion- like. For the phonon mode, the peak shape is 
determined by the function B{e/U,e' /U). It has a maxi- 
mum dX e' = U (e < U). In diffusion-like modes there is 
a divergence (integrable) at e' ~ e (see (|104|l and (lllll) ). 
These general forms of the curves are the same for stain- 
less steel and Fomblin, but the relative contribution from 
phonon-like and diffusion-like modes are quite different for 
the two substances. Contributions from phonons and diffu- 
sion into small energy transfer can be estimated from (|92|l 
and H104() . (|lll|l . Omitting numerical factors, one has for 
diffusion-to-phonon contribution ratio 



Wd 

w. 



p 



(118) 



For thermo-diffusion, the high value of the ratio c/vq 
may be compensated by the smallness of the quantity 
a = Cp/C'v ~ 1 and additionally suppressed by the factor 
l/V^D, when do is large. This is just the case for stain- 
less steel, when the right part of (|118|) is near unity. As for 
Fomblin, the right-hand side of (|118|l is near 40. This qual- 
itative estimate agrees with numerical results presented on 
Figs. m 13 and 13 The thermo-diffusion contribution is com- 
parable with the phonon one for stainless steel (Fig.^ and 
dominates for Fomblin fiuid (Figs. |21 and inj. 

The transverse relaxation mode in liquids for small en- 
ergy transfer is also of diffusion type, which originates from 
sound waves with high damping due to viscosity. The con- 
tribution of this mode, as compared to the longitudinal 
phonon one, can be evaluated from (|118|l with a = 1. It 
far exceeds unity in agreement with Figs. |21 and |31 



It should be emphasized that the two diffusive modes, 
which, due to the large ratio H118|l . give the main contribu- 
tion to the small energy transfer for Fomblin, are governed 
by different physical parameters. The transverse mode has 
no small parameter a, but its quantity d d^, Ij 109(1 dif- 
fers from that for thermo-diffusion, do (|103|l . by magnitude 
and temperature dependence (see Table 1). 

Remark. The transverse mode was considered here with 
a simple liquid model ((35|I - H37|I . where transition to "ideal" 
liquid with very low viscosity is not correct. In this paper 
we have performed calculations for d > 1 (see Sec. lA. 3.111 . 
so the definition (|109|l gives the lower bound for the vis- 
cosity ly. (In fl^ we have argued, that for small d < 1 
the contribution of diffusion mode to inelastic scattering 
vanishes.) 

Let us now discuss quantitative results for the proba- 
bility per bounce Hll()|) for the neutron to change its en- 
ergy from U up to (10 — 10"^) • U, that is to go from 
ultra cold to very cold energy region ("small heating"). 
For stainless steel, the probability for neutron transition 
to the region from U to (10 — 10'^) • C/ is of the scale of 
~ 10~^^ — 10^^^ (Fig.©, while for Fomblin it is of the scale 
of ~ 10-^° - 10"^ (Figs. Eland El). Evidently, these contri- 
butions could not be seen in the measurements |SlEllZllHlin|- 
Nevertheless note, that our calculations for Fomblin shows 
a definit increase of low-energy heating with the tempera- 
ture due to the contribution from diffusion-like transverse 
mode. It is in a qualitative agreement with the data |B1IS|- 
The measured effect of small heating on Fomblin seems to 
be explained by surface effects (see jHUHl). 

As a general conclusion, phonons are ineffective for small 
energy transfer both for stainless steel and Fomblin. Ther- 
mo-diffusion is even less effective for stainless steel, but 
gives very large effect (one-two orders of magnitude larger 
than phonons) for Fomblin, where the transverse diffu- 
sion mode gives comparable contribution and even exceeds 
thermo-diffusion when the temperature goes up. (Note, 
however, that it goes down more steeply when the trans- 
ferred energy becomes higher.) 

In the region of high energy transfer, the ratio of thermo- 
diffusion and phonon (at G = 0) contributions can be es- 
timated from and H106|l . 



Wd 

w. 



ac 



p 



vod 



D 



(119) 



It is of the scale of 10^^ for stainless steel and 10^ for 
Fomblin fluid. Thus, thermo-diffusion may dominate over 
phonons not only for small energy transfer, but as well for 
a large one, provided the dimensionless thermo-diffusion 
coefficient d^i is low. This conclusion may have sense only 
for liquids and amorphous materials, since for crystals, and 
stainless steel in particular, the main process for the tran- 
sition up to the temperature region, upscattering, is that 
with momentum exchange with the lattice, G ^ (see 
SecEiSl- 

Therefore, the processes responsible for high energy 
transfer are quite different for crystals and amorphous ma- 
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terials. Our result for coherent UCN upscattering on pure 
crystalline material H99() coincides with that previously 
obtained. It is known that its estimate (|117|l is one-two 
orders of magnitude too small to be observed in a real stor- 
age experiments. The reason is, evidently, a much higher 
absorption rate due to the radiative capture and possible 
upscattering on hydrogen contamination of the surface. 

For Fomblin only the surface effects have been discussed 
in details so far. We have not found any published re- 
sults with evaluation of escape probability from the bulk 
properties of this material. In other words, there is no un- 
derstanding why the parameter / (see ^01) of bulk losses 
lies in the range 10^^ to 10~^. Thus our result just for bulk 
effect given by the integral (|ll(i|l (calculated with Sf = £_d) 
is new. Numerical results for Fomblin presented in Figs.El 
and El give for the loss probability per bounce (~ /) the 
following values: 0.5- lO"'^ at 293 K and 0.7- lO""^ at 373 K. 
Possible reasons for the discrepancy of one-two orders of 
magnitude are mentioned in conclusion. 

9 Conclusion and outlook 

A theory for coherent inelastic scattering of ultracold neu- 
trons was developed with a quite general form of the 
correlation function, which allows to consider energy- 
momentum exchange of the neutron with sound waves 
(phonons) as well as with diffusion- like relaxations. The 
latter appears when a correlation function is considered in 
the hydrodynamic region of momentum and energy trans- 
fers (the low limit for q and lo). Namely, we discussed 
the longitudinal thcrmo-diffusion mode both for solids and 
liquids as well as the diffusion-like transverse mode for 
liquids. 

In the frame of this theory, we evaluated the probability 
per bounce for an ultracold neutron to gain energy and 
to go into cold or thermal region. In this paper the case 
G = (i.e. no momentum transfer between neutron and 
the crystal lattice) is consistently considered for the first 
time. 

Numerical estimates are presented for the stainless steel 
storage bottle and one with Fomblin-fluid coated walls. 
For Fomblin, the most explored experimentally, theory 
predicts the overwhelming effect for thermo-diffusion pro- 
cesses, two-three orders of magnitude larger than from 
phonons, and a quite strong dependence on the temper- 
ature originated from viscosity. This result is in a qualita- 
tive agreement with experiments. 

However, the numerical values for the upscattering prob- 
ability per bounce obtained for stainless steel as well as for 
Fomblin are still one-two order of magnitude too small to 
fully explain experimental data. Thus, it is worthwhile to 
mention restrictions of the presented results and possible 
future improvements. 

Results from diffusion processes are quite sensitive to 
the parameters of correlation function. It is very desir- 
able to clarify the structure and parameters of correlation 



function, especially for liquids. In fact, it can be done by 
experiments in optics and other means outside the neutron 
physics. 

We considered our samples as a spatially uniform media. 
This may not be correct for Fomblin, which exhibits some 
polymer features and may have a non-uniform structure at 
the scales comparable with the ultacold neutron reduced 
wave length ~ 10 nm. To get the correlation function for 
polymer and to include it in the neutron scattering theory 
is a serious problem for the future. Nevertheless, the effect 
may be evaluated from simple physical arguments. The 
probability for energy transfer, as was mentioned above, is 
proportional to the intrusion length or, what is the same, to 
the time of neutron-wall collision. Spatial polymer struc- 
ture changes, in particular, the self-diffusion parameters 
for the neutron and increases the time, spent by the neu- 
tron inside the sample during the collision. This time can 
be evaluated from transmission experiments and the corre- 
sponding correction can be easily introduced into the for- 
mulae of this paper. 

On the other hand, gigantic polymer molecules can be- 
have as nano-particles. If it is the case, elastic scattering of 
the neutron by the whole moving molecule in their center- 
of-mass system would result in inelastic scattering of neu- 
tron in the laboratory system. Similar mechanism of small 
cooling and heating due to nano-particles at the surface of 
solid materials was proposed in )3Ul . 

No surface effect was included, which seems to be of 
no importance in solids for UCN large wave length. This 
argument may not concern specific surface excitations in 
liquids, say, capillary waves 22 , 10 , which can be included 
into consideration as one of the excitation modes. This 
particular case requires a different correlation function and 
is the subject for a separate paper. 

Several important effects were not included in this pa- 
per. Among them are spin-flip effects and surface con- 
tamination with hydrogen etc. Note, for example, that in 
0] evidences are presented that the total upscattering of 
UCN is mostly related just to hydrogen contamination of 
the walls of material traps. All these effects are of impor- 
tance for incoherent scattering, not included in this paper, 
which deals only with coherent inelastic scattering of ini- 
tially ultracold neutrons. 
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A Calculation of inelastic cross sec- are real and positive. The function fl{q±) takes the form 

tion and inelastic transition prob- _ i 

^^.j.^y niq^) = ^^iUiq^,z*)-Uiq^,z)), (125) 

The aim of this Appendix is to simpUfy the general expres- ^^^^^ z = x + ty and 

sion (|62(l " (|63|l for inelastic cross section. The first stage is ^ ^* 

calculation of the integral ESI. U{a, z) = H ; — . (126) 

a — z* a + z 

A.l Calculation of the integral Jit{\v) This complex function of two complex variables wiU be 

' extensively used below and it is worthwhile to exhibit its 



symmetry 

C7(a, z) = U{-a, -z) = -U{-a, z*) = U*{a*,z*). (127) 



We start with a general integral of the form 

/(A ri) = — / dq± A(gx — A)A(gj^ — ri)!F{q±) (120) ^he real and imaginary parts of this function for the first 
J argument P + are of the form 

— oo 

2 fix 

where A and rj are complex parameters, and ^{q±) is a Re C/(/3 + ij, z) — + (7 + '^uY ~ kP) 1 (128) 



function with poles above ((7„|) and below {qni) the real 



axis, but not near it. The function Af.xl is defined bv 16011 . , ^ 2x , o , s , , ^^^^ 

andtherefore ImU{(3 + ^j, z) ^ -- {^^j + {j + 2y){x' + {j + yf)) , 

(129) 

. , , , . , , 1 where 

A(,, - A)A(,, - ,) = X A = + + (7 + yry ~ 4/3^x^ (130) 

/ (A - 7/)a i(q^~^)a -i(q^-^)a\ '^ote, that U{ij, z) is pure imaginary and given by 

X 2 cos eV ^/-eV - 

We transform (|120|l into the contour integral closing a 

path in the upper or the lower half-plane. Then, neglecting Now, we are ready to apply the general result (|122{ and 

the term proportional to e-i™«"Ta and e-|i™«"il'', we get P23|l to the integral It is enough to find three "ele- 

for any positions of A and 77 mentary" integrals /„(A, rj) (n = 1, 2, 3) with 

2 sin(A- 77)0/2 J'i{q±)^n{q^), 
/(A, r;) = (.F(A) + Hv)) 

^ (122) ^(9±) 

2cos(A-??)a/2 ^2{q±) =2,2' (132) 

+ ^ -— (ViX) ~ V(j])) , 91+91 ^ ^ 

A — 77 

e function ") 

of.F(q±) 



g±»(g±) 
•?^3(9±) = — Y- 

where the function V{X) is defined by the residues in poles + gj^ 



All J>i(qj_) have the poles from f2(g^), the last two have 

( , r\ the additional poles in = izgii. The functions Pri(A) 

y res^ _ y res^ ^^^^ ^^^^ 

To apply the resuhs for (jHSJ, one needs to specify the '^^^^^ " ~4|Ip ^(•^'^)) ' 

pole structure of the integrand. For the factors Fi^t{q±), 

it is evident from ^ and The function f^i,t(g±) _ ^i(A) _ »A^(zg||,z) 

is given by H41|l and can be presented as the sum of four ' q'^ + A^ 2(7|||zp(g^ + A^)' 

simple pole terms. The poles are located symmetrically " " 

with respect to the axes at the points x ±iy and —x ± 

where ^ A^ 2|z|2(q2 + a^) ' 



^/p^Tf^ - p2 / y^p4 _^ p4 _^ p2 Note, that J^n and P„, which determine /„(A, 77) by (fT^ . 



(124) 



contain the function U{a, z) p26|l varied by first argument. 
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The function J(A, 77) (|65|l can be presented as a sum of 
J„(A,r?) with coefficients evident from and ljH7|l . After 
some rearrangement the results can be written in the form 



(134) 



Both terms in the right-hand side are proportional to the 
factors U {a, (3), however, J^^^ includes all terms with a = X 
or r], while j'^' includes terms with a = zgn . Namely, 



. {X — r})a . (X—r])a 



X-rj 



Ji.ti^,v)^cos $i,t(A,?7) 



where 



A/,t(A,?7) 



911 1 



FiAX)U{X,z)~Fu{f])U{r],z* 



(135) 
(136) 



(137) 



*i.t(A,?7)=± 



(qi + X^)(qj + rf) 



?f(Q||q||+Agr )(Q||q[|+?7Qr ) 



- (AQiiqii - QT gpi^Qiiqii - QY if)^ 

(138) 

and the functions Fi^t are given by l|6H|l and 1)67(1 . In ((138|l 
"+" and "-" correspond to "f and respectively. 

Note, that the dependence on the parameter a, the thick- 
ness of the layer, is explicitly indicated in (|135|) and p36|) . 
All other functions are independent on a. The functions F 
and $ are constructed entirely from kinematic factors, ini- 
tial and final neutron momenta. Parameters of correlation 
function, z = x + iy, apart from the explicit factor 
enter only in the function U{a, z) with a = A, 77 or iq^^. 

The cross section (|62|l . (|63|) is the function of the in- 
coming and outcoming neutron wave vectors, which are 
included in parameters Q|| = k|| — kj| (|52|l and Q^i^ , QY 
((54|l . For the crystal these parameters are split into the 
parts inside the first Brillouin zone, q|| , A, 77, and corre- 
sponding reciprocal lattice vectors Gii, Gx, G^. Note, that 



Fi{X) = Fiir,) = Qj + QY QY 



0, ^uiX,ri)^±Ql 



(139) 
(140) 



Ft{X)=Ft{7^) 

when all G||, Ga, G,, are zero. 

Thus, with calculated J{X,r]), the general expression for 
the inelastic cross section is of the form (I62|) . (I63II as the 
sum over the indices ct, tr', r, t'. This summation, for a 
general case, is cumbersome because the variables A and 77 
((64|l enter the functions U{X, z) and U{ri, z) in a sophisti- 
cated way. 



A. 2 Inelastic cross section for initial sub- 
barrier neutrons 

In this part of Appendix we specify the general expres- 
sion for the inelastic cross section (|62f) . H63|) for the case, 
when the initial neutron is sub-barrier. The final energy 
can be both below and above the barrier. The scheme of 
approximation varies for these two cases, and they will be 
considered separately. The object of analysis will be the 
factor Ili^t in the cross section (|62ll . All other factors are 
universal and will be added afterwards. 

For the initial sub-barrier neutron it is useful to define 



k — z3e, 86 



4-2 i.2 



(141) 



The quantity 7 = g-*-" g for the initial channel may 
be neglected everywhere. In particular, for the neutron 
coming from the left we have 



0. 



(142) 



To compensate 7 in the first factor in H63|) we need from 



the sum 7 



gffia_ Therefore, the factor needed comes 



only from the first exponent in H135I) with a = t = 1. The 
same exponent should be taken into account in ((136|l . 

Following ((134|l . the 11 factor H63|l can be presented as 
the sum of two terms 



n 



(1) 

i,t 



X (y)-'-72(y*)-r72 



X-r, ■ 



(143) 



n 



lniU{iqii,z) ^ 



(144) 



X (7')-'^'/2(y*)-r72 $^V(A,,7). 



Parameters A and 77 are defined by ((64|l and can be rewrit- 
ten now as 



A = ias -1- a'k[, 77 



(145) 



where k[ = k' — G± is the part of k' in the first Brillouin 
zone (G_L = cr'Gx = r'G^). The identity e*'^^" = 1 for 
the reciprocal lattice vector is taken into account. The 
functions A and $ in 1(143(1 and ((144(1 depend on indices a' 
and t' via their arguments. 

A. 2.1 Final state: below the barrier 

To be more precise, only normal to the surface component 
of the final neutron energy should be below the barrier. 
For the normal component of the final neutron wave vector 
inside the sample we introduce (with analogy to 1141() ') 



k' 



88 



1,2 y2 



(146) 
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and 



In the sums (|143f) . (|144|l over a' and r' the main term On the other hand 
corresponds evidently to cr' = t' = 1 because the others ^ ^ ^ ^ 

contain the factor 7' = e"*'", which may be neglected ^i,t ^^i^*) = ^i,t (~'^0' ~-^o) = A;_i(Ao, Ap), (154) 
apart from a small vicinity of the barrier, where k'^ —^ka, 
se'a < 1 and 7' may approach the unity. The barrier region 
is discussed in details later. Note that even for a' = t' = 1 A,^j^^(A, —A) = A^^ ^^^(Ao, — Aq), 
only the coefhcients A''-^'' contribute to (|143|l . H144|l . but 

not A'^^\ Surely, it means that the neutrons coming from A^^^^ ^(A, —A) = A^^^^ Aq, A^) = —^]^^ *(Ao, — Aq). 
the left are scattered also to the left. ' ' ' (155) 

For the small energy transfer one may put G,, = and ^i^j^ ^he help of ^^^^^ one readily finds 
Ga = G,, = 0. Taking into account (|139|l and p4U|) . we 

get 1 / ImA"(Ao,A^) 

A,"(Ao,A5) ^ (Qf + |Aonim^(Ao,z) ' \l - r'WY \ 2a. ^^^^^ 



2A„ 2|A„||zP ' +2Re A',,A'^iy)-^^h±ShlZ^ 

n«=o, V 2Ao 

(147) __ 

Im[7(iQ||,z) Q|| Imf/(iQ||,z) The sum m a similar way results m 



n,, -$,,(Ao,Ao) 2Q |,|2 -± 2\z 



(148) nf,)^ ^"^^(^g|h^) ((l + /^)$ii(Ao,Ag) + 
where 29,1 jz^l - r'^'^p 



Xo =i{iB + «;'). (149) 



+2Re (A;iA'!i(7')-^<i>,V'(Ao,-Ao))) 



-2 
(158) 



Note, that for pure imaginary first argument a = 17 the (157) 

function ImU is given by ^here the quantity <^j\{Xo, X^) CM is real. 

The last terms in ()i56|l and (|157|l have different values 

A. 2. 2 Final state: above the barrier for backward (B) and forward (F) scattering because 

In the sums H143|l . (|144|l over a' and t' all terms should be r^, . /n-it s _ / _ // i\ 

considered equally important, since 7' = e**^ " is now not a L +1 -iVM J +1 -1 w 7 Kf ) 

small parameter. If we neglect radiative capture, then the i(R) /(r,)* 

quantities fc' and r' are real, and 1 7' I = 1. [^^1^-1(7') ] = ^+1 (7') = -r' ■ (159) 

It is useful to split H143() , (|144|l into the sums with a' = t' 

(where ?? = A*) and that with a' = -t' (where 77 = -A) A. 2. 3 Final state: near the barrier. Resonances 

of penetration 

^ / A''^'''^'fA A*) 

n^"'^-' = I X^M' ,p — - + With the results of two previous subsections we are able to 

^■^ 11 /2 /2|2 \ / J ' cr ' n - 

|i ~ 7 I \ cr' 2*86 consider the upscattering of UCN into the regions close to 

jy<j' (^x — A)\ barrier, below or above. 
+ '^^A'^, A'*^,{Y)~'^ ] , In the low vicinity of the barrier ae' ^ ae. Above the bar- 
er' 2A J j-jg], "vicinities" should be distinguished, "wide vicin- 

(150) ity", where k' < k^, but k'a > 1, and "close vicinity", 
lmU{iq\\,z) f , 2 a'cr' * where fc'a < 1. 

^i,t ~ 2q bPll — r'^7'^P i -^^^"^'^ ^^^^ In the "wide vicinity" of the barrier, ae' and k' may be 

" ^ treated as small parameters. Due to small energy transfer, 

A' A'* ( '\-<^'(b''' (\ -\\\ ^^"^^ Ga = G^ = 0. In t he a pprox imati on A ~ Aq = «ae 

2^^T'^-<y'y^ I ^i,t y-^' j ■ and 77 Ag the equations l{T15|l and ((Till can be written 

(151) All 



For the explicit summation in H150|l we need to know IlJ ^ ~ ^{r',j') — '■ — , (160) 



(1) ^ „^^, Ai,t(Ao, Ag) 

2«as 



how A;^j'"'(A, A*) and A;^|-'''(A, -A) depend on a'. From 

(|137|l with the use of symmetry properties (|127ll one can ^ ^j-^,' 'y')<j)ii(;\(j X*'j lrn£/(zQ||, z) 

show that the quantity ' ' 2Q|||z 



S(/,y)<&i".(Ao,A5) ^^VJ'J ^ (161) 
Ai](A,A*)^Aij(Ao,A5) (152) -here the factor 

is pure imaginary, where -(r',7') = ,, _ ,2 ,212 X! ("^'^ ''^^(7'*) 

' ^ ' (t',t' 

Ao = ia3 + A:;. (153) (162) 
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for backward and forward scattering is equal to 



= (r , 7 ) 



|2 _ ^n^a\2 ' 



^ (r ,7 ) 



(163) 

When the final momentum k' is just above the barrier, 
where it is reasonable to put r' ~ 1, then 1 and 

0, and the corresponding barrier values of Hi^t ()160|l 
and (|161|l coincide with that from (|147|l and H148() . The 
deviations arise only close to the points where 



y2 ^ g2ifc'a 



1 



k'a — Trn, 



1,2,3. 



and, therefore, a difference between r' and 1 is of impor- 
tance. For these energies 11^ decreases and 11^ increases 
up to the peak values 



n 



(1)B,_F 
l.t 



1 A,",(Ao,Ag) 
4 2isi 



„(2)B,F_1 11,. Im[/(iQ||,z) 



(165) 



(166) 



It is seen from (| 164(1 that these minima in and maxima 
in arise when the thickness of the layer is divisible 
by the half of the neutron wave length in the substance. 
Thus, the reason for the oscillation of the cross section is 
the resonances of penetration through the layer. 

When the final momentum is not too close to the bar- 
rier, the factors containing 7' — e**^ " (e.g. |l — r'^7'^| ) 
strongly oscillate. With increase of the neutron final energy 
the amplitude of oscillations goes down because r' 0. 
However, in the "wide vicinity" of the barrier the oscil- 
lations are of importance. For not very thin sample these 
oscillations cannot be resolved and only the smoothed cross 
section is of interest. 

A. 2. 4 Smoothing over oscillations 

The functions Ft 1(156(1 and 1(157(1 with account for ((158(1 
and l(159() oscillate due to the factors 



^ = V 



^/2(/+m) ^2i(l — m)k a 



(167) 



^/2(/+m) ^2i(l — m — l)k' a 



Lm=0 



Keeping only the terms with zero indices in the exponents, 
we get the smoothed functions 



1-r'- 



ImAi](Ao,Aa) 
2a: 



2Ao 



(168) 



(n 



Lt I 



Im U (ig[| , z) 
2g||kP(l-r'4) 



(l + r'2)$ii(Ao,AS) 



2r'"Re$;\-i(Ao,-Ao) 



(169) 

where n = 3 for backward and n — \ for forward scattering, 
and Ao is given by 1(153(1 . 

The smoothing is not justified in the "close vicinity" 
of the barrier where k'a ~ 1. But since the smoothed 
quantities 1(168(1 and 1(169(1 have the same limiting values as 
1(160(1 and 1(161(1 . one may use the smoothed cross section 
everywhere above the barrier. 

The functions A/_t, which enter ((168(1 . are disclosed from 
inU, (ESJ and as 

ImAii(Ao,AS) = 



Im 



■ (Qiiqii + Ql'Ao)(Q||q|| + Q^l *Ao) C/(Ao, z) 



kP(9f + A§) 



(170) 



Im Aii(Ao, AS) = - Im 1^ \q\ + \Q^l\^ - 

(Q||q||+Qi^Ao)(Q||q||+Qii*Ao) \ [/(Ao,z) ' 



Ai-i(Ao,-Ao) 



((Qiiq||)'-(gl^Ao)2)[/(Ao,^) 



(171) 



kP(9f + A§) 



(172) 



Ai-i(Ao,-Ao) = - I Q\~{Q'tf- 

(Q||q||)'-(Ql'Ao)n C/(Ao,z) 



9^ + Ag 



with 



Q^^ — i£e + k' , Ao = iae -I- k[. 



(173) 
(174) 

The functions which enter ((169(1 . are evident from 

(tim . 

A. 2. 5 Interim conclusion 

The general formula for the inelastic cross section ((62(1 for 
the initial sub-barrier neutrons is now fully defined with 

Ui^t = ^\^t + U[^t given by (fTT7| and (dHl for the final 
state below the barrier {k'j_ < kg) or by 11168(1 and p69(l for 
that above the barrier {k'j_ > ko). 

A. 3 Inelastic transition probability for 

G = 

In this part of Appendix we calculate the dimensionless 
factor ((70(1 , which determines the inelastic transition prob- 
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Figure 7: Integration area for the double integral H175(l : 
the cases k' < fc|| (a) and k' > fc|| (b). 



ability H69() with momentum transfer inside the first Bril- 
louin zone (G = 0). 

In the integral (|70|l it is useful to replace the variable k'j^ 
by fcjp and the variable (/9 by Q| = (ky — ky)^ running from 
(fcJi - fc||)2 up to (fc[| + Thus, 

WiAk±,k\\^e') = 

k^ f f (175) 

= /dfcf /dg^5(Q||,fc[|,fc||)i^nz,t, 



where the function g can be presented in two forms 

1 



5((5||,fcj|,fc||) = 



^((fc||+fc||)2-Q|) (gj|-(fc||-fc||)2 



Note, that 



U'|2 _ / 4fc^/fco, _ 



4fc^V(fc^ + fc')', k'^>ko. 



(176) 



(177) 



Integration area for double integral (|175|l is shown in Fig.[7| 

The general formulae for the factors 11;. < are given by 
(|147|l . (|148|l for the final state below the barrier (fc^ < ko), 
and by (|168|l . H169(l for the final state above the barrier 
{k'j_ > ko). When G = 0, the Eas. ((TC5|) and are sim- 

plified, and the sums of forward and backward scattering 
contributions can be written in the form 



n 



(2) 
Lt 



± 



Q|| ImC/(i(5||, ; 
(l + r')|z|2 



(179) 



n 



(1) 



1 



(Q2 + |Ao|2)Im{/(Ao,z) 



1 - r' 



£C\Z\ 



(178) 



where Xq = iee + k' . 



A. 3.1 Estimations for variables and functions 

The function U{a, z) H12t)|l in (|178|l and (|179|l represents 
the correlation function. The real and imagine parts of the 
variable z ~ x + iy are determined by the equations 



(180) 



With the fixed initial and final neutron energies, e and 
e', the "damping factor" (see SectionOJ may be treated 
as the constant (defined by appropriate physics) . The same 
is true for the quantity f.uP' jt? . It is convenient to use one 
of these parameters as a unit and transform all variables 
to the dimensionless form (k!^^ fc||, Qy Qy, flj^t II;_t 

etc.). We choose as the unit r/V2, thus, e.g., 



X —* X 



\/2a 



V2y . ^/2Q|| 

-, c^ii ^ y4\\ - 



and, in particular. 



r 



r2c2 



r 



(181) 



(182) 



Note, that e = 1, except for the longitudinal thermo- 
diffusion mode (see text after (|41() '). where e = 0. In the 
integral H175|l the upper limit for fc'p is 



2t'2 

r2 ■ 

The equations (|18()(l take the form 

1 ~2 ~2 /S2 -2 

xy = \, y - X ^ Q» - uj . 



(183) 



(184) 



They define the functions a;(Qji) and y{Q\\) provided the 
parameter cj^ jg given. Therefore, for our problem (with 
fixed £ and e'), there is only one additional parameter, a)^, 
which should be specified by the physics of the process. 

The damping factors for longitudinal sound waves 
both for solids and liquids (|24l) and for transverse sound 
waves for solids (I25II . as a rule, are much less than k'^. On 
the other hand, in the diffusion-like cases (|34|l and H37() the 
damping factor takes the form 



p2 ^ H ^ 2m|e' 
D 



2mD 



(185) 



where D = Ds for longitudinal thermo-diffusion both for 
solids and liquids, and D — v for damping transverse waves 
in liquids. Parameter d is dimensionless. It means, that 
1183|l can be presented as 



2e'd 



(186) 
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Note, that in all realistic cases d ^ 1 (see Section [T)). The factor (|147|l . H178|) for the small energy transfer 

Therefore, both for sound-like and diffusion-like modes one takes the form 

~„ , , nP^-S- (198) 

> 1. (187) ' 

The following relationship But it cannot be substantially simplified for large energy 

^ transfer. However, this is of no importance. Indeed, in this 

^ ^ ! (188) f,a,se the main contribution to the cross section comes from 

, „ , . , T 1 1 the region (I) due to the dominating term 119511 . 

is also of great importance. Indeed, o w o ' 

^ ^ _£ £ yj2 ^ fl89l A. 3. 2 Estimation of the integral 

Let us now estimate the integral (|175|l . We see that the 

because e' is restricted from above by the Debye energy region (I) contributes only to the factor Wi for the large 

2 \i/3 energy transfer due to H195|l . Since A:' '--^ Qy ae ~ , one 

— c{ IT n) , [ ) perform the integration over fcy in 1)175(1 in the narrow 

and Ed < 2mc'^ for all substances. interval Q|| - < < Q|| + fc|| (see Fig.|7|D). It gives 
Let us separate two regions of integration oveq Qy in 

in3: (I), where Q fc'» landed (see CH3 and CHi), w['\k^, k« ^ e') ^ ^ /dg^^Hp), (199) 

therefore, 27rfc J H fc^ 
(I) f<l<y-Q||-fc', (191) 

and (II), where Q,, « k' , therefore (see CHZJ and UHHJ), ^^""^^ ^ ~ % is given by Taking 

into account ((184|l and ((191|l and integrating over (3| from 

(II) i,y,Q||«fc'. (192) ^2 to ^'2 nne obtains 



In the region (I) the factor 11; ' ((148(1 and p79(l is of the ofc 1 T^fc'^ 

n!J~ji~3i«i- (193) ' 

II where the result is presented in terms of the physical quan- 

The explicit form of the factor ftp'' in the same region titles. 

depends on whether small or large energy transfers. For Now let us consider the contribution to the cross section 

small the energy transfer, y ~ Q,, ~ fc' ~ ae, and one has l™i^ ^he region (II). It is convenient to rewrite the factor 



« 1 (194) 



1(175(1 as follows 



^\\ M^,,(fc^,fc|l_.£') = ^ JdQl^Jdkf~gf{k'^), 

both from p47(l and 1(178(1 . The case of the large energy (201) 

transfer corresponds to the relations: ae ^ fco <C y ^ Q|| ~ where the function Ili^t/(f>' given by ((196(1 and 1(198(1 de- 

k'. In this case only the region fc^ > k^ is of interest, pends only on Q|, and 
Therefore, one obtains from 1(178(1 

/ N <P'\t'\^ 4 / , 

~ 4 (qi + k'A fik'^) ^ = 71 fcl - (^ik'l - kD^k'i - ki 

~ ^. (195) k'^ kl\ ^ 

se (k'l + yA (202) 

with e{x) = for X < and e{x) = 1 for x > 0. 

Note, that this term dominates due to the relatively small The integral over fc'p is of the form 
factor ae in its denominator. 

Now let us consider the region (II). The dimensionless /" ~ „ - 4 / _ fz ^ \ 

r . fV(2)u.^,u, J u A u • u v dk\?gf(k',) = ^{l(k')-I(Jk'^-k^)] , (203) 

factor / both below and above the barrier may be writ- J \\ y ^' 1^2 y ^ ' v/j^ \ / 

ten in the form 



jj(2)^^^;!Wl±m^, (196) 



where 



I^P (S' + {Q\\+ y?) ' I{b) = 0{b^ - fc|) j dkf ~g e{b^ - fcf)^62_fc/2^ (204) 

J fc^ < fco, The variable fc||^ runs from (Q|| — fc||)^ up to (Qii -I- fc||)^ (see 

I {k'j_ + k')/k'j_, k'^ > fco- ^^^'^'^ Fig. [7J|. Then, I{b) can be presented in terms of elliptic 
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integral 



7r/2 





(205) 

where 



a2=4Q||fc||, /32 = 62 _ (g^i _ 



(206) 



It is of importance, that in the region (II) both for small 
and large energy transfer one has 



fell < fc'fc|| < - fc^, 



■ll^-l 



a < /3. (207) 



Then, really E{a'^/P'^) ~ 7r/2, and, therefore, the contri- 
bution from the region (II) to the cross section is 



^.l\2 



(208) 



where k', 



fc'2 - fcjj, while n, 



(1) 



and nj ^ are given by H198|l and l|196|l , respectively. 

To calculate the integral it is convenient to replace 
by X with the use of 



A. 4 Inelastic transition probability for 

G ^ 

In this part of Appendix we calculate the dimensionless 
factor (|73|) , which determines the inelastic transition prob- 
ability (|72|l with the momentum transfer G to the crystal 
lattice. The final neutron energy e' ~ fi^G^ /2m is much 
larger than the barrier energy U and initial energy e, thus 

k|| ~ -G|| , ~ fcl ~ G_L > 33, r' -> 0, i' -> 1. 

(215) 

The n factors (|168() and (|169() have now equal contribu- 
tions from backward and forward scattering and, with the 
help of approximations H215(l and expression Ij 170(1 - (|174|l . 
H138|l . take the form 



n 



(1) 



1 / (G||q|| +GxAo) ?7(Ao,0) 
Im 



<ll + >^l 
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1 ^ / (GiiQii +GiAo) t/(Ao,z) 
Im 



(216) 
(217) 



-Zf + Ag 



dQ 



(209) 



The limit Q\\ — corresponds to 



(210) 



In the opposite limit Qy — > oo, i — s- 0, then the integral 
((208|l can be presented in the form 



W^"\k^,k^^^e') = !^f{k'^)(j^'^-J^'^ 



Wrik^,k\\^e') = ^.nk'^)j('^ 

TTk 



where 



+2y) 



dx. 



(211) 
(212) 

(213) 



„(2) _ ^ gf|G||q||+G^Aop-|AoG||q||-G^g|| 



where 



kf + Agp 

Im U{iq\\ , z) 
Ao=ia3 + ^_L, ^^ = k'^-G^. 



(218) 
(219) 



The integral (|73|l is over the Brillouin zone with the cen- 
ter in G. It is convenient to take 



d-'k' = dk'^d^k\ d£,^d^q\\, 



(220) 



because qy = — k'| — G||. Integration over the angle be- 
tween qii and Gii gives 



2„2 



(G||q||)=0, ((G[|q[| 



Then the W factors take the form 



(221) 



Straightforward calculation gives 
J 



2) _ 2xo 



(214) 



Then, summing the contributions from the regions (I) 
and (II) we get in terms of the physical quantities the re- 
sults (TT ^ -fTS Ii . 



W^f{k^,k\\^G) 



d(,±dqu X 



q^ {Gj /2-Gi 



TTkxG^ J "'^"""ll 

^U(Xo,z) 



Im 



9^ + Ag 



G^^ II ImC/(Ao,z) 



(222) 
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G) = ±- 



TTkG^ 



d£_±dq?: X 



911 



(q/2-Gi) (g||-|Aop)Im£/fa||,z) 



(223) 



Now we can integrate over ^± transforming the inte- 
grals into contour ones by closing a path in the complex 
plane. For three functions to be integrated, one easily 
obtains 



V In + \ 



0, 

TT Im U{iq\\ , z)/q\\ 



j Im U{Xo, z) = ~2ttx, 



0, q\\ > se. 



9|| < ae, 
9|| > ac, 
(224) 

(225) 
(226) 



Using these results we get for the sum of and 
factors 

Wi,t(fci, fc|| ^ G) = T ^ ^ ^ X 



X / dq 



9|| Im;7(zg||,z) 



(227) 



2fc| 



ka^G^ J "|z|2 
Using (|131|l and (|209|l . one obtains 

dgfll^'",^^""^=-4j(^), 



(228) 



= 2JW, 



where J^^^ and J^^^ are given by 1214|) . Thus, the W fac- 
tors are of the form 1)79(1 and ((8()|l . 
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